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First VLBI mapping of circumstellar ^^SiO maser emission 
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Abstract. We report the first VLBI map of the v=0 7=1-0 maser line of ^'SiO in the long-period variable star IRC + 10011. 
We have found that this maser emission is composed of multiple spots distributed in an incomplete ring, suggesting that this 
maser is also amplified tangentially, as already proposed in other SiO circumstellar masers. We present also VLBI maps for 
the 7 mm ^*SiO v=l and 2 7=1-0 and the 3 mm v=l 7=2-1 lines. The ^'SiO masing region appears to be located in a layer in 
between the ^**SiO v=l 7=1-0 and ^**SiO v=l 7=2-1 lines. In addition, we confirm that the 86 GHz maser v=l 7=2-1 forms 
in an outer region of the circumstellar envelope compared to the other ^^SiO masers studied. Finally, we discuss the possible 
implications of the observational results on the SiO maser pumping theory. 
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1. Introduction 

SiO maser emission is found in the innermost shells of the 
circumstellar envelopes of Long-Period Variable (LPV) stars, 
evolving along the Asymptotic Giant Branch. To date, three 
different isotopomers of the SiO molecule are known to ex- 
hibit maser emission towards these objects: ^^SiO, ^'^SiO and 
^°SiO, with relative abundances of [2^SiO]/[2''SiO]~ 20 and 
[2'*SiO]/[^"SiO]~ 30 respectively. 

Strong ^^SiO maser lines have been detected, either 
by single-dish or interferometry techniques, in hundreds of 
evolved stars, from the v=0 to the v=4 levels and up to the 
rotational 7=8-7 transition (e. g. Pardo et al. ,1998) . In con- 
trast, the ^^SiO and ^°SiO circumstellar emission has been less 
studied. In fact, maser amplification from these isotopomers 
has been measured in a few sources and a few rotational tran- 
sitions. The first detection of a ^'SiO maser line, the v=0 
7=1-0, in a variable star was performed by Cho et al. ('1986i. 
Subsequent single-dish studies of the v>0 7=1-0 ^"^SiO cir- 
cumstellar masers revealed that this rare isotopomer emission 
also had some properties usually associated to the ^^SiO lines, 
such as the correlation with the IR 8 /urn radiation or its time 



variability (Alcolea & Buiarrabal ll992l l. 

Our theoretical knowledge of these rare isotopomer emis- 
sions is still poor. Some pumping mechanisms have been pro- 
posed, but they do not compare well with observations. In order 
to constrain these models, we have studied several SiO maser 
lines in the OH/IR variable IRC +10011 (WX Psc) by means 



of Very Long Baseline Interferometry (VLBI) techniques. We 
present in this paper interferometric maps for the 7 mm v=l, 
y=2 7=1-0 and 3 mm v=l 7=2-1 transitions of ^"^SiO, as well 
as the v=0 7=1-0 line of the ^'SiO isotopomer. For the first 
time, we have been able to detect and map the ^''SiO maser line 
at high spatial resolution (better than 1 mas). We focus on the 
comparison between the maps of the different maser transitions 
and discuss how these results may affect the overall theoretical 
SiO pumping scenarios. 



2. Observations and data analysis 

Using the NRAO' Very Long BaseUne Array (VLB A) we per- 
formed sub-milliarcsecond resolution observations of the SiO 
maser emission in the O-rich LPV star IRC +100 11 on 2002 
December 7. 

The results in this paper correspond to the second epoch 
of observations for IRC +10011, as part of a VLBA multi- 
epoch/transitional study of the SiO maser emission in AGB 
stars (see for further details Soria-Ruiz et al. i2004j . The ob- 
served maser lines are summarized in Table Q] In this obser- 
vational run all the SiO transitions were detected and mapped, 
including the ^^SiO line. This is the first time that maser emis- 
sion other than the main molecular species has been imaged 
using very long baseline interferometry. 

The data were correlated in Socorro (New Mexico). In the 
43 GHz observations, a bandwidth of 8 MHz was recorded and 
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Table 1. Summary of observational results 
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Fig. 1. The total power (continuous lines) and recovered 
flux (dotted lines) spectra of the SiO masers observed in 
IRC+10011. The dashed line indicates the systemic velocity 
of the source (Vlsr^ 9.2 kms"', taken from Cernicharo et al. 
[19971 . 



the correlator provided 256 frequency channels, thus achieving 
a spectral resolution of about 0.22 km s ' . For the 86 GHz tran- 
sition, the bandwidth was 16 MHz and 512 spectral channels 
were used, being the resolution ~ 0. 1 1 km s"^ The calibration 
was done using the standard procedures for spectral line obser- 
vations in the Astronomical Image Processing System (AIPS) 
package. For the final phase calibration we have used a strong 
maser component as reference, and therefore, all positions in 
the resulting maps are refered to the location of this spot, which 
is not the same for all the transitions. 

Fig. n presents for each of the SiO transitions, the total 
power spectrum of Los Alamos antenna, which was used as 
reference in the data reduction process, and the spectrum of the 
flux recovered in the maps. The Une profiles of the v=l and 
v-2 7=1-0 masers are similar although they diff'er from the 
other two transitions. About 75% of the v=l and v=2 7=1-0 
maser emission was recovered after the calibration and imag- 
ing, while 10% and 25% of the emission was imaged for the 
^'^SiO v=l 7=2-1 and ^'^SiO v=0 7=1-0 respectively. 

The integrated intensity maps are shown in Fig. |2] As can 
be seen in the upper panels, the spatial distributions of the 
^^SiO v=l and v=2 7=1-0 emissions appear to be similai", both 
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are ring-like, with the v=2 located in a slightly inner layer of 
the envelope. In contrast (see lower panels of Fig.|2}, the ^^SiO 
v=l 7=2-1 and ^^SiO v=0 7=1-0 maser maps are composed of 
a fewer number of spots and the ring distribution is less clear In 
order to study and compare more accurately the angular extent 
of the different masing regions, we fitted a ring to our observa- 
tional data. To do this, we selected only the maser features with 
a signal-to-noise ratio larger than 6 in at least three consecutive 
spectral channels. These calculations give the center, the mean 
radius (R) and the width of the ring (aR, being twice the stan- 
dard deviation of the sample). We have summarized in Tabled 
for the transitions observed, the inner and outer radii derived 
from the fits (defined as Ri^ - R - jAR and Roui = R + \aR) as 
well as the shape of the gaussian restoring beam at half power 
(major and minor axis and position angle of the major axis) of 
the resulting maps. 

3. Results 

Previous long baseline interferometric studies of the SiO 
masers demonstrated the existence of very compact emission 
distributed in ring-like structures (see e.g. Diamond et al. 119941 
and in IRC +1001 1 by Desmurs et al. l2000> . Those results are 
limited to the ^^SiO transitions. Therefore, our new detection 
and mapping of the ^'SiO v'=0 7=1-0 illustrates for the first 
time the spatial distribution of this rare isotopomer emission in 
evolved stars. 

The map of the ^'SiO transition is composed of seven maser 
spots with velocities ranging from 5.7 to 11.3 kms The 
brightest feature has an intensity of ~1.77 Jy beam^'-kms"' 
(see bottom-right panel of Fig. 2). This ^^SiO emission is the 
weakest among the four lines observed, with an intensity 6 
or 7 times lower than the ^'^SiO 7=1-0 transitions. The emis- 
sion forms a ring, though incomplete, with a mean radius of 
-13.5 mas, therefore this maser radiation is probably amplified 
tangentially, as other well studied ^'^SiO maser transitions; in 
Table [Owe show the inner and outer radii obtained from our 
maps following the method described in Sect. 2. 

Both the v=l and v=2 7=1-0 ^**SiO emissions present a 
ring-like structure composed of 10-12 spots with a similar dis- 
tribution. In IRC +1001 1, Desmurs et al. f^U^ and Soria-Ruiz 
et al. ( 2004) found a systematic shift between these two emit- 
ting regions of about 1-3 mas, with the v=2 radius being al- 
ways smaller. We confirm this result since the v=2 ring ra- 
dius is also ~1 mas smaller than the v=l one (see Table [2 and 
Fig. |2j. We note that a similar trend has also been observed 
in other Mira variables in several epochs over a stellar period 
(Cotton et al. .2004. Yi et al. .2005.1 . From the fliree existing 
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Fig.2. Integrated intensity maps for the SiO maser emissions in IRC+10011. Top; ^ SiO v=l (left) and v-2 (right) 7=1-0 
rotational transitions. The peak flux is 9. 14 and 9.72 Jy beam ' ■ km s ' respectively, the contour levels in both images are 0.28, 
0.57, 1.14, 2.28 4.16 and 8.34 Jy beam"' -kms"' and the rms noise level 0.07 Jy beam"' -kms"'. Bottom: ^^^SiO v=l 7=2-1 
(left) and ^^SiO v=0 7=1-0 (right) maser transitions. The peak flux is 3.37 and 1.77 Jy beam"' • kms"' respectively. The contour 
levels are 0.15, 0.20, 0.40, 0.80, 1.60 and 3.20 Jy beam"' -kms"' for the 7=2-1 line, and 0.08, 0.12, 0.23, 0.46, 0.92 and 1.76 Jy 
beam '-kms"' for the 7=1-0 line. The rms noise levels in the maps are 0.02 and 0.01 Jy beam"' kms"' respectively. Circles 
denote the ring fitting of the masing regions (dashed: mean radius R, continuous: Ro^t and defined as ^ + \^R)- We note that 
the four images have the same scale to ease the comparison between the emitting regions and that in all the cases the derived 
center of the ring is placed in the center of the box, while the origin of the map remains at the position of the reference spot (see 
Sect. 2). 



maps of these 7 mm lines in IRC +1001 1 (Desmurs et al. 120001 
Soria-Ruiz et al. 120041 and this work, which correspond to 
phases, ~ 0.8, 0.1 and 0.5 respectively) changes of about 10% 
and 20% are derived for the sizes of the v=l and v=2 7=1-0 
maser emitting regions respectively. 

From our maps of rotational transitions within the same 
vibrational state, fliat is, flie ^**SiO v=l 7=1-0 and 7=2-1, 



we conclude that the latter maser is produced in a shell fur- 
ther away than the 7=1-0 one (see Table 0. Furthermore, the 
spatial distributions of the various components clearly differ 
(Fig-EJ- A. similar relative location of the ^'^SiO v=l 7=1-0 
and 7=2-1 masing regions has been found in previous obser- 
vations of IRC +1001 1 (Soria-Ruiz et al. 2004). These authors 
have pointed out that there is a clear contradiction between the 
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observational results, that is, the 3 mm transition in the v- 1 
being produced in an outer region of the envelope than the cor- 
responding 7 mm one, and the theoretical predictions, in which 
these two lines are spatially coincident (either radiative or colli- 
sional models). The authors conclude that this discrepancy can 
be overcome when the overlap between ro-vibrational transi- 
tions of H2O and ^^SiO are introduced in the ^'^SiO excitation 
(see Soria-Ruiz et al.|2004i. Other comparisons of these 3 and 
7 mm lines have been done by Doeleman et al. (2004 1 in the 
Orion KL nebula, finding a similar result, and by Phillips et al. 
J2003> towards the AGB star R Cas. In this case, the authors 
claim a similar distribution of both emissions. 

The ^^SiO v=0 7=1-0 masing region appears in a layer lo- 
cated in between those of the ^^SiO v=l 7=1-0 and the ^'^SiO 
v=l 7=2-1 emitting shells. This is in fact remarkable since 
these two y=l transitions require excitation temperatures of 
about 2000 K and the ^'^SiO maser is a low-excitation line 
[£(7=1) ~2 K]. Therefore, the obtained location of this mas- 
ing region suggests that the rotational levels involved in the 
inversion are populated mainly via de-excitations from upper 
v>0 states. 

4. Discussion 

Attempts to model the '^SiO maser amplification in evolved 
stars, and in particular the v=0 7=1-0 line, have been made 
using different excitation schemes. Robinson & Van Blerkom 
<1981j and Deguchi & Nguyen-Quang-Rieu ( 1983j proposed 
that the ground state ^''SiO masers are produced if the vibra- 
tional transitions of this molecule present a significantly higher 
opacity along the radial direction than in the tangential one. 
This asymmetry with the direction is produced in slowly accel- 
erated envelopes, yielding line profiles composed of two peaks 
separated by twice the expansion velocity of the envelope, sim- 
ilar to what is typically observed in circumstellar OH masers. 
These predictions are clearly in contradiction with the obser- 
vations, since the line shapes of ^^SiO masers are composed of 
narrow peaks near or at the stellar velocity (see also Alcolea 
& Bujarrabal 1992 1. Moreover, the map for IRC -i-lOOl 1 shows 
a ring-like geometry, very probably indicating tangential am- 
pUfication (as it happens in the ^^SiO maser Unes). Finally, a 
small velocity gradient is not expected for the inner region of 
the envelope where the ^^SiO maser is produced. 

The other proposed mechanism involves line overlaps be- 
tween infrared transitions of ^^SiO and ^^SiO. This effect was 
first suggested by Olofsson et al. C1981.) to explain the v'=0 
7=2-1 ^^SiO maser. Subsequent calculations by Gonzalez- 
Alfonso & Cernichai-o (1997"t suggest that the v=0 7=1-0 
^''SiO line can be efliciently pumped by the overlap between the 
^'^SiO v=2-l 7=4-3 and the ^'^SiO v=l-0 7=1-0 ro-vibrational 
transitions. 

These models do not provide calculations of the size and 
brightness distribution of the emitting regions, thus making it 
very difficult to conclude on the compatibility of the theory 
with the presented VLBI observations. 

In any case, our observational results are in agreement 
with the predictions of the mechanism proposed by Gonzalez- 
Alfonso & Cernicharo (.1997; . For this model to work, the 



ro-vibrational transitions v-2-l of SiO must have relatively 
high opacities (t> 1). This condition is also required for the 
arising of the strong ^'^SiO v-1 and 2 masers (in both colli- 
sional and radiative excitation schemes). Therefore, the fact 
that the -^SiO emission arises in the same inner shell as the 
^^SiO v=l and 2 masers guarantees that the main requirement 
of the model by Gonzalez- Alfonso & Cernicharo (1997 1 is sat- 
isfied. (We also note that the ^''SiO v=l-0 7=1-0 line also coin- 
cides in frequency, better than 2 km s"' , with the ro-vibrational 
line V2-2-1 72,6-83,5 of H2O.) Finally, as it has been men- 
tioned in the previous section, overlaps also appear to be re- 
sponsible for the discrepancies found between observations of 
some rotational lines of ^'^SiO and models (see also Bujarrabal 
et al. 1996 Herpin & Baudry 2000 1. 

We are currently carrying out similar high resolution stud- 
ies of maser lines of ^**SiO and ^'^SiO in other evolved stars, to 
see if our results are reproduced in other targets, and especially 
the observed location of the emitting region and spatial distri- 
bution of the 29siO v=0 7=1-0 line. A more extended sam- 
ple will also contribute to a better understanding of which pro- 
cesses may be involved in the pumping and if the line overlaps 
play such an important role in the SiO excitation. 
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